Deeply inelastic pions in the exclusive reaction p(e, e'ir + )n above the resonance region 
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A model for the p(e,e'-ir + )n reaction which combines an improved treatment of gauge invariant 
meson-exchange currents and hard deep-inelastic scattering (DIS) of virtual photons off nucleons 
is proposed. It is shown that DIS dominates and explains the transverse response at moderate and 
high photon virtualities Q 2 whereas the longitudinal response is dominated by hadronic degrees 
of freedom and the pion electromagnetic form factor. This leads to a combined description of the 
longitudinal and transverse components of the cross section in a wide range of photon virtuality Q 2 
and momentum transfer to the target t and solves the longstanding problem of the observed large 
transverse cross sections. The latter are shown to be sensitive to the intrinsic transverse momentum 
distribution of partons. 
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At Jefferson Laboratory (JLAB) the exclusive reaction 
p(e, e'ir + )n has been investigated for a range of photon 



virtualities up to Q l 



5 GeV at an invariant mass 



of the Tt + n system around the onset of deep-inelastic 
regime, W ~ 2 GeV 0, 0, ||. A separation of the cross 
section into the transverse <tt and longitudinal ctl com- 
ponents has been performed. The longitudinal cross sec- 
tion ctl is well understood in terms of the pion quasi- 
elastic knockout mechanism Q because of the pion pole 
at low —t . This makes it possible to study the charge 
form factor of the pion at momentum transfer much big- 
ger than in the scattering of pions from atomic elec- 
trons @. On the other hand, the ctt is predicted to be 
suppressed by ~ 1 /Q 2 with respect to ctl for sufficiently 
high Q 2 »A^ CD Q. 

However, the data from the n-CT experiment [|| show 
that cr T is large at JLAB energies. At Q 2 = 3.91 GeV 2 
ctt is by about a factor of two larger than ctl and at 
Q 2 = 2.15 GeV 2 it has same size as ctl in agreement 
with previous JLAB measurements [l|. Theoretically, 
the model of Ref. [7|, which is generally considered to 
be a guideline for the experimental analysis and ex- 
traction of the pion form factor, underestimates ctt at 
Q 2 = 2.15 GeV 2 and at Q 2 = 3.91 GeV 2 by about one 
order of magnitude [3J] . Previous measurements at values 
ofQ 2 = 1.6 (2.45) GeV 2 Q show a similar problem in the 
understanding of ctt- Even at smaller JLAB 2] and much 
higher Cornell @ values of Q 2 there is a disagreement be- 
tween model calculations based on the hadron exchange 
scenario and experimental data; see Ref. [t| for a possible 
interpretation and references therein. 

In this work we first generalize the treatment of Ref. Q 
for the longitudinal contribution. We then propose a res- 
olution of the ctt problem. The idea followed here is to 
complement the soft hadron-like interaction types shown 
in Figure [T] which dominate in photoproduction and low 
Q 2 electroproduction by direct hard interaction of vir- 



tual photons with partons followed by the hadronization 
process into n + n channel, to form the 7r + - electropro- 
duction framework. As we shall show, then the large ctt 
in the reaction p(e,e'7r + )n can be readily explained and 
both ctl and ctt can be described from low up to high 
values of Q 2 . 

The exclusive reaction 



e(P e ) + N(p) - e'{P' e ) + Tr(fe') + N(p') 



(1) 



with unpolarized electrons is described by four structure 
functions ctt, ctl, ctlt and cttt [ 101 ] - After the integra- 
tion over the azimuthal angle between the leptonic and 
hadronic scattering planes only ctt and ctl remain and 
the differential cross section takes the form 



da e / 1 dQ 2 dudt 



7T$ 



E e (E e -u) 



[daT/dt + eda L /dt},(2) 
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where e is the virtual photon polarization. The defini- 
tion of the virtual photon flux $ follows the convention 
of Ref. [l(J. The subscripts T and L denote the projec- 
tions of the (7* , 7r) amplitude onto the basis vectors of 
the circular polarization of the virtual photon quantized 
along its three momentum q: T- transverse (A = ±1) and 
L— longitudinal (A = 0) polarizations. 

At first we consider the soft hadron-exchange part of 
the 7r + -electroproduction amplitude. In Figure [T] the 
Feynman diagrams describing the high energy tt + electro- 
production in the hadron-exchange approach are shown. 
It has been well known for a long time that the 7r-pole 
amplitude, first diagram in Figure [T] gives the dominant 
contribution to the longitudinal response ctl. The tt- 
pole amplitude by itself is not gauge invariant and charge 
conservation requires an addition of the electric part of 
the s-channel nucleon Born term (third diagram in Fig- 
ure [T]) 0, [Tl| • When considering realistic vertex func- 
tions, which include form factors, current conservation 
is violated and one has to restore the gauge invariance 
of the model [l2[. A simple solution to this problem is 
to choose all of the electromagnetic form factors to be 
the same 0- However, it is experimentally known that 



FIG. 1: The diagrams describing the hadronic part of the it - electroproduction amplitude at high energies. See text for the 
details. 



these form factors are not the same and have different 
scaling behavior -F 77r7r ~ 1/Q 2 for the pion form factor 
and Ff ~ l/Q 4 for the proton Dirac form factor. 

In the following we use the Regge pole model of Ref. Q 
which is based on the same set of Born diagrams but 
concerning the electromagnetic form factors we employ a 
prescription proposed in Refs. [HI, [3] where an arbitrary 
form factor F(Q 2 ) can be accommodated by the following 
replacement of the currents 

-» r' M (Q 2 ) = r» + [f(q 2 ) - 1] v?r u , (3) 

where = g^ v — q^q v jq 2 stands for the projector into 
the 3-dimensional transverse subspace. This procedure 
guarantees that the resulting current r lfi obeys the same 
Ward-Takahashi identities as _T M . Thus, and as long as 
gauge invariance is implemented for real photons, one 
can use the experimentally determined form factors in 
the 7r-pole J% and s-channel nucleon Born currents 
and still retain gauge invariance for arbitrary Q 2 . 

Making use of Eq. ([3]) the 77T7T and jNN vertex func- 
tions are given by 

rj^ = e(k + k'Y + e [F 1WW (Q 2 ) - 1] V^(k + k') v , (4) 
^nn = e^ + e [jf (Q a ) 1] V^lu, (5) 

where the four momentum vectors of pions are k (incom- 
ing) and k' (outgoing). In Eq. Q we have - as usual - 
assumed that the half-off-shell form factor F 77r7r (Q 2 , i) 
depends only on Q 2 . Using Eqs. (Tj| and (J5J) the gauge 
invariant hadronic current J M describing the reaction 
p(7*, n + )n is constructed as a sum J M = J% + J£. 

At high energies the exchange of high-spin and high- 
mass particles has to be taken into account. To account 
for these states we replace in J£ the 7r-Feynman prop- 
agator by the Regge propagator Furthermore, since 
the s-channel Born term can generate the pion pole it- 
self we factor out the pion propagator in the sum 
J% + Jj^ following Ref. [7[ and reggeize it according to the 
above prescription. The hadronic current which satisfies 
the current conservation, i.e. q p J^ = 0, takes the form 



-iJ^ = V^g^NN u s /(p')75 



+Ff(Q 2 ) 



(k + k'y i 

t-ml + i0+ 
W 2 - Ml + i0+ 



x[t-ml + i0 + ] 
where 



+ [f-MQ 2 )-f!{Q 2 )] 

2a»(i) 



(k - k'Y 



Q 2 



Us(p) 



w 
mi 



-i-KOL^ (t) 



sin(7ra 7r (i)) T(l + a w (t)) 



al + a'j = 0.7(t - ml ) 



(6) 



(7) 



is the degenerate 7r-6i-trajectory, Wq = 1 GeV and the 
Gamma function T suppresses the singularities in the 
physical region (t < 0). In Eq. (J6]) g^NN — 13.4 is the 
pseudoscalar irN coupling constant, t = k 2 , k = kl — q = 
p~p' and other notations are obvious. It should be noted 
that in the current ^ the two different form factors of 
the nucleon and the pion appear; this is in contrast to 
the work of ,7f where these two form factors F 17r7I and 
Ff were assumed to be identical in order to reach gauge 
invariance. For the pion charge form factor we use a 
monopole parameterization 



F 1 ^{Q 2 ) = {l + Q 2 /k 2 y 



(8) 



with the cut-off A 77r7r as a fit parameter. The Dirac form 
factor Ff(Q 2 ) is described by a standard dipole form. 

The second diagram in Figure [T] describes the exchange 
of the p-meson Regge trajectory. The current J£ reads 



iJ» = -iV2G pNN G 1P KF 1P AQ 2 )E^ af:i qvk 



xu s >(p') 



W 
Wo 



(i + Kpb/3 

2a p (t)-2 , 



2M, 



iP + p')p 



-i-KClp(t) 



sin(7ra p (t)) T(a p (t)) 



.00 



(9) 



The parameters needed for the proper description of the 



current J£ are 



a p (t) = 0.55 + 0.8* 



(10) 



as the degenerate p-a2-trajectory, G p nn — 3.1 is the 
vector and n p = 6.1 is the tensor pN coupling constants. 
The 7/57T coupling constant G 7p7r = 0.728 GeV -1 has 
been deduced from the decay width [1 51 ] 



+7T" t "7 



67.5 keV. 



(11) 



For the ■jmt vertex form factor F lpv we use the prediction 
of Ref. fijl. 
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FIG. 2: The longitudinal dai^/dt (top panels) and transverse dtJT/dt (bottom panels) differential cross sections of the reaction 
p(j*,ir + )n at average values of Q 2 = 1.60 (2.45) GeV 2 [l[ and Q 2 = 2.15 (3.91) GeV 2 0]. The solid curves are the contribution 
of the hadron-exchange model and the histograms are the contribution of the DIS pions. The discontinuities in the curves 
result from the different values of Q 2 and W for the various — t bins. The dashed histogram in the lower left panel shows the 
contribution of the DIS pions for the average transverse momentum of partons yj {k 2 ) = 0.4 GeV. 



In Figure [2] the results for the p("f* ,ir + )n differen- 
tial cross sections da^/dt (top panels) and da^/dt (bot- 
tom panels) are compared with the data from JLAB 
Ftt-2 Q and n-CT [H experiments. The longitudi- 
nal cross section da-^/dt is very well described by the 
hadron-exchange model (solid curves) with the cut-off 
in the pion form factor being fixed to the constant value 
A 2 ,^ = 0.52 GeV 2 . The discontinuities in the curves 
result from the different values of Q 2 and W for the vari- 
ous —t bins. The steep fall of da^/dt away from forward 
angles comes entirely from the rapidly decreasing 7r-pole 
amplitude. The interference of this amplitude with the 
s-channel nucleon Born term is minimized due to the 
presence of different form factors for both amplitudes. 
The contribution of the natural parity p-exchange is neg- 
ligible in ctl and or- 

This comparison with data shows that F 77rw can indeed 
be reliably extracted from the longitudinal data. 

Again, the model strongly underestimates da^/dt, for 
example, at Q 2 = 1.6 GeV 2 by a factor of 10 and at Q 2 = 
3.91 GeV 2 by a factor of 30. This is also seen in the model 
of Ref . , although somewhat less pronounced [H, Q . 

The solution to this problem is still missing. One might 
describe this transverse strength in the language of per- 
turbative QCD by considering higher twist corrections to 



a GPD based handbag diagram. However, such a calcu- 
lation does not exist and it is not clear if a higher-twist 
expansion converges in the kinematical regime consid- 
ered here. Our solution of this problem, therefore, is 
to model such effects. We start from the observation 
that the second term of Eq. ^ contributes very little 
to both the longitudinal and transverse cross sections. 
Here only the nucleon Born term is taken into account 
to conserve the charge of the system. However, at the 
invariant masses reached in the experiment (W ~ 2.2 
GeV) nucleon resonances can contribute to the l7r chan- 
nel. Similar to the replacement above of the pion propa- 
gator by a Regge propagator that takes higher meson ex- 
citations into account we now complement the s-channel 
nucleon Born term with direct hard interaction of vir- 
tual photons with partons (DIS) since DIS involves all 
possible transitions of the nucleon from its ground state 
to any excited state [17|. Note, that our suggestion con- 
cerning the partonic contribution follows the qualitative 
arguments in [l8| where it has been shown that the typ- 
ical exclusive photoproduction mechanisms involving a 
peripheral quark-antiquark pair in the proton wave func- 
tion, the i-channel meson-exchange processes considered 
above, should be unimportant in the transverse response 
already around Q 2 > 1 GeV 2 and play no role in the 
true deep inelastic region. This we have already seen in 
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Figure H 

The total transverse DIS cross section reads 

DIS _ ATT 2 aFP(x,Q 2 ) _ ATT 2 aFP(x,Q 2 ) Q 2 + 4M p 2 i 2 
° T ~ 1-x uM p ~ 1-x Q 4 l + K(x,Q 2 )' 

(12) 

where a ~ and x — „H f . In the following we as- 

137 2vM p & 

sume that a partonic description of deep-inelastic struc- 
ture functions F^ 2 ^{x,Q 2 ) works well not only in the 

Bjorken limit where 1Z = a^ ls /a^ ls tends to zero but is 
valid down to values of Q 2 considered in Figure O 

To determine the structure of events in DIS a model 
for the hadronization process is needed. Furthermore, 
since at JLAB (Bjorken x > 0.3) the antiquark con- 
tent of the structure functions becomes negligible, we 
model the DIS by the j*q — > q knockout reaction fol- 
lowed by hadronization through string fragmentation. In 
the present description of hadronization in DIS we rely on 
the Lund model (LM) [19| as depicted in Figure [3] where 
the 7*q — > q process followed by the fragmentation of an 
excited colored string (wavy curve connecting the quark 
lines 1 ) into two particles (irN) is shown. The LM pre- 
dicts two jets for the 7r + n final state in the forward and 
backward directions. As a realization of the LM in DIS 
we use the pythia/jetset implementation [2(|. The LM 
involves parameters which have been tuned in different 
fragmentation channels. Our approach here is to modify 
as few parameters as possible compared to the default set 
of values [20J which describes the ir + SIDIS spectra mea- 
sured at JLAB 2l| remarkably well. Since in pythia the 
average transverse momentum of partons y (&t ) cannot 
be fixed from first principles and since it affects the slope 
and magnitude of dar/dt at forward angles we choose 
this as a free parameter. Therefore, one has to regard 
the average \J (k 2 ) used here as an effective parameter 
which is tuned to obtain an agreement with data. How- 
ever, for consistency we use the same value for y (k 2 ) in 
all kinematic regimes together with the default jetset 
parameters. As pointed out above we view the string 
fragmentation process as an effective model for higher- 
order twist effects, for example in GPD based handbag 
calculations. The success of our description may then be 
taken as an indication that the string fragmentation pro- 
cess described in jetset works well down to the rather 
low invariant mass of about 2 GeV where the individual 
nucleon resonances tend to disappear. 

The lower part of Figure [5] shows that da^/dt receives 
the dominant contribution from DIS fragmentation pions 
(solid histograms). In Eq. (Tl2l) for F% we use the fit of 
Ref. [12] and for 1Z the parameterization of Ref. [23| has 
been employed. In [2(j the value of \J (fc 2 ) = 1.2 GeV has 
been used for all Q 2 bins; this value is close to the default 
pythia value of \J (k 2 ) — 1 GeV and well within the com- 




FIG. 3: A schematic representation of the partonic part of the 
7T- electroproduction mechanism. The wavy line represents a 
color string. See text for the details. 

mon range of transverse momentum distributions [24[. 
As one can see in Figure [2] (bottom panels) the absolute 
value and the —t dependence of da^jdi are very well re- 
produced. A decrease of y (k 2 ) increases the slope and 
magnitude of d<JT / dt at forward angles. In Figure [2] this 
is shown for \J (fc 2 ) = 0.4 GeV (dashed histogram). 

We have also compared the model results with data 
from the JLAB Ftt-1 experiment Q at lower values of 
Q 2 . Also here we find that an addition of the DIS pions 
describes the experimental data very well. However, con- 
trary to the situation at higher values of Q 2 where the 
hadronic part gives only a marginal contribution to ctt, 
at low Q 2 the problem of double counting arises when 
using both the DIS and the Regge contributions to the 
transverse cross section. Following Ref. [25[ this could 
be solved by turning off the leading order DIS contribu- 
tion, as required by gauge invariance for "f*q — ► q, when 
approaching the photon point where the Regge descrip- 
tion alone gives a good description of data Q- In the 
calculation presented here the transverse part is solely 
generated by the DIS process (TT2"]) without any further 
modification. 

In Figure 2] we confront the result of our calculations 
(solid curves) with the new JLAB data [26| for unsepa- 
rated cross sections at average value of W — 2.2 GeV. 
The data are very well described by the present model 
in a measured range from Q 2 ~ 1 GeV 2 up to 5 GeV 2 . 
Furthermore, assuming that the exclusive cross section 
behaves as a^ IS {Q 2 ) cx Ff(x, Q 2 ) in Eq. JT2J) 2 and that 
the ratio 1Z is small or nearly the same both for protons 
and neutrons we predict then a smaller transverse cross 
section in the reaction n(e,e'n^)p off neutrons, i.e. 

a^/a p T ~ Ff/ff « < 1, (13) 

while because of the rr-pole dominance 

<l< * 1. (14) 



1 Not to be confused with the perturbative one gluon exchange. 



2 This behavior has been already noticed in @] from a fit to data 
and is supported by the present model. 
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FIG. 4: Differential cross section da/dt — d<JT/dt + edai^/dt of the reaction p(7* 
predictions. The experimental data are from Ref. [261 ] . 



The solid curves are the model 



A preliminary analysis in Ref. [l| has shown that the 
latter ratio is indeed consistent with unity and ctl/ct 
must be larger for ir~ than for tt + . This, together with 
the fact that ctl is described very well also at the highest 
Q 2 by the Regge picture alone indicates that the DIS 
contribution to the exclusive longitudinal channel must 
be small. 

In summary, in this work we have extended the ear- 
lier model of Ref. [7j such that the electromagnetic form 
factor of the pion and the nucleon no longer have to be 
set equal in order to achieve gauge invariancc. In addi- 
tion, we have proposed a resolution of the ctt problem 
in the reaction p(e, e lr K + )n above the resonance region. 
A model which combines the gauge invariant hadron- 
exchange currents and DIS of virtual photons off partons 
has been proposed. The model with hadronic states as 
the active degrees of freedom describes the longitudinal 
cross section ctl very well and exhibits the dominance 
of the 7r-pole mechanism while ctt is grossly underesti- 
mated. We have shown that the description of cjt at 
values of Q 2 > 1 GeV 2 requires a proper inclusion of the 
hard scattering processes and that 7*5 — > q followed by 



the TT + n fragmentation of the nucleon may naturally ex- 
plain the large transverse cross section observed at JLAB. 
The model can be used for the extraction of the pion 
form factor from high energy pion electroproduction data 
with longitudinally polarized photons. The sensitivity 
of the transverse cross section to the transverse \J (k 2 ) 
of partons can be used to reduce the theoretical uncer- 
tainties in the interpretation of the color transparency 
signal observed at JLAB in the reaction (e, eV + ) off nu- 
clei [27 , 28] . Finally, we mention another ctl /ctt puzzle in 
the reaction p(e, e' K + )A(T0) (29j which may apparently 
get a similar solution [301 ] . 
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